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Abstract 

We study the influence of the space and time dispersion on the fre- 
quency dependence of the wave vectors of electromagnetic waves propa- 
gating in three-dimensional photonic crystals. Two types of structures are 
considered: media with weak periodic modulation of the permittivity, and 
photonic crystals composed of the periodically arranged identical resonant 
dielectric particles. It is shown that in these systems, in contrast to elec- 
trons in solid crystals, different types of excitations exist. For example, 
a peculiar kind of polaritons arises in the photonic crystals due to the 
interaction of the electromagnetic field, eigenoscillations of the dielectric 
medium, and Debye resonance. The widths of the transparency zones and 
of the band gaps have been calculated as functions of the frequency and 
of the parameters of the media. It is shown that in the photonic crystals 
with dispersion, the number of transparency bands is larger than in non- 
dispersive systems, and the width of the gaps in the frequency spectrum of 
photons depends on the wave vector. The interaction of different types of 
waves deforms the Brillouin zone, so that it may not have a plane boundary 
(for example, a sphere), in which case the classical Bragg condition does 
not hold. 

PACS: 78.20.-e, 71.36,+c, 71.35.-y 

1 Introduction 

The propagation of waves in periodic structures is a topic of profound impor- 
tance, both for understanding the general properties of matter-radiation interac- 
tions, and for the ever-increasing number of practical applications (for example, 



1 



photonic crystals) [1], [2], [3]. In spite of the close similarity to the quantum- 
mechanical problem of electrons in a crystal lattice, the propagation of classical 
waves in photonic crystals exhibits many rather unusual distinctive features. The 
most dramatic example of such a difference is the effect of time and space dis- 
persion, which has no analogue in the motion of conduction electrons. While the 
propagation of electromagnetic waves in homogeneous dispersive media has been 
thoroughly studied [4], the effect of dispersion on the radiation propagating in 
periodic structures still remains poorly understood. The exceptions are periodic 
multiple quantum well structures (periodic arrays of dielectric layers), for which 
the polariton spectrum and the optical response have been studied in considerable 
detail [5], [6], [7], [8]. In the recent publication [9], the exciton-polariton suscep- 
tibility matrix was calculated, taking into account the well-well interaction, and 
the polariton dispersion curves for an infinite periodic array of quantum wells 
have been computed. Some interesting polariton-induced effects were discov- 
ered. However, the one- dimensionality of the system considered prevented the 
investigation of important effects inherent only in dispersive media of higher di- 
mensionality. In the present paper, we deal with the linear wave equation of the 
most general type and study the spectral properties of dispersive periodic dielec- 
tric media, namely, the effect of the time and space dispersion on the frequency 
dependence of the wave vector. The methods employed (small perturbation and 
local perturbation methods) impose no restrictions on the dimensionality of the 
system. As formulated, the problem is closely related to the calculation of the 
energy of an electron as a function of the wave vector in periodic crystals. How- 
ever, in contrast to the later case, the solution of the corresponding dispersion 
equation for electromagnetic waves can have several branches, i.e. different types 
of waves can propagate. The interaction of these waves is the crucial factor in 
the formation of the photonic spectrum, and gives rise to rather unusual effects. 
In particular, it deforms the Brillouin zone, so that it may not have a plane 
boundary (for example, a sphere), in which case the classical Bragg condition is 
violated. 

We consider the most general form of the wave equation 

^idf' idi) ^ idr 1 idt^ ^' ^ ^ 

d d 

where H and h are arbitrary functions of — - and and h is a periodic 

ior iot 

function of r with periods r n = riifi + n 2 r 2 + n 3 f 3 , where fi : r 2 , r 3 define the unit 
cell of the photonic crystal, and ni, n 2 , n 3 are integers. In non-dispersive media, 
Eq. (1) becomes a classical second-order wave equation; space-time dispersion 
corresponds to higher-order derivatives with respect to the coordinates and the 
time in Eq. (1). 

According to the Bloch theorem, the solution of this equation has the form 
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E(r,t) = (2) 

where E{f) is a periodic function with periods r n . 

Our goal is to investigate the spectrum of waves propagating in the periodic 
system, i.e. to find the dependences, u>i(k), of the frequency ui on the wave vector 
k (the subscript / indicates different types of waves). The periodic functions E{r) 
and h(f, J^) can be expanded in Fourier series, 

E(r) = E n e 1 ^, (3) 

n 

^^-) = E e " V M^,-), (4) 

n 

— * 

where k n - are the reciprocal lattice vectors multiplied by 2%. 

Substituting Eq. (2) and (3) in Eq. (1) and setting the coefficients of each 
e t>Cmr equal to zero, we obtain the following infinite system of equations, 

H m (x m ,uj)E m + h m - n (x„,uj) E n = (5) 

n 

—> — * 

with x m = k + k m 

This system of equations has a nontrivial solutions, E n , if its determinant 
vanishes, 



D(x m ,u) = (6) 

— * 

Equation (6) is the dispersion equation that determines the function cu(k). 
From Eq. (2), it follows that u(k) is a periodic function of k with periods k n 
[10], [11]. To avoid problems associated with the infinite rank of the determinant 
in Eq. (6), we consider two limiting cases: a weak periodic modulation that 
corresponds to the nearly-free electron approximation, and a crystal, for which 
the method of local perturbation is applicable (an analog of the tight-binding 
approximation). 



2 Weak periodic modulation 



In this section, we assume that the amplitude of the periodic modulation of the 
medium is small, so that 



and the term h(r, ^JU, -^)E in Eq. (1) can be considered as a small perturbation. 

When h = (zero-order approximation) all coefficients /i m _„ in Eq. (5) are 
zero, E m = 0, except for m = (£o = const), and the dispersion equation for 
the homogeneous medium takes the form 



H(k,<j)=Q (8) 



Note that only real solutions uii(k) correspond to propagating waves. We 
denote the number of such waves by L. To find coefficients E n to first order of 



the small 



h 



, there is no need to deal with the exact dispersion equation, Eq. 

(6). Instead, we set to zero all E n with n ^ in the second term in Eq. (5) and 
obtain 

E n = jz (9) 

H [k n , uj 

(10) 

It follows from Eq. (7) that \E m \ « \E \ except for the case when there 
exists some m = v for which H{k v , w{k v )) is small. In this instance E v in Eq. (5) 
is not small compare to E , and one has to keep in Eq. (5) also terms containing 
E v . This yields 

H(k,u?jE Q + h- v {K_ v ,oj)E v = (11) 
H (k,u) Eu + h v (x u ,u)E = (12) 

The condition for nontrivial solutions of Eqs. (11), (12) to exist leads to 
the following dispersion equation in the periodically modulated medium in the 

degenerate case 



H (jt, u?j H (xj,, cu) — h v (k, uj^j h~ v (£- v , oS) = 0. (13) 
To solve the dispersion equation Eq. (13) we write H (k, u?jm the form 

L 

H(k,u) =f(k,u) H[^-^i(k)] , (14) 



=i 



where function / has no poles on the real axis. Substituting Eq. (14) into Eq. 
(13), we obtain 

l J_ hJk,u)h„Jk,u) 
l\ L - cu t (k) J] [co - u m (k v )\ v — )_ = (15) 

1=1 m=l / (k,uj f{Z v ,u) 

If we write the solutions of Eq. (15) as 

u (jcj = oji (Pj + 8ui (Pj , (1 = 1, ....L) 
and assume that for all non-interacting waves 

uji (Pj - Ui (kjj » 8ui (Pj , 



(16) 



(17) 



the corrections 5u>i ( k ) can be found by successive approximations. This proce- 
dure yields: 





h v (k v ,ui (Pj 


) 


to 


f (k, ui (Pj 






Ui (Pj - U rn (Pj 
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- u m (kjj 



l = 



(18) 

It is important to note that this correction is of second order in the small 
perturbation h. 

The situation is different when for some I = n, there exists k for which u n (k) 
is close to uj n (k u ), and inequality Eq. (17) is violated. Then, the following 
approximate equation for the dispersion curves uj n {k) can be derived from Eq. 
(15): 



uu — u) n [k 



[oj - u n (h v )\ -r?Ak,K v ) = 



(19) 





h v (k v ,u> n (k 


) 


2 


/„ (k, u n (Pj 




- (%,)] nf=i 


UJ n - LOi (k 





The two solutions of the quadratic equation Eq. (19) are: 



(20) 



uj n [k v + u> n [k 
u n± {k) = W ± 



u n [k) -co n (k v ) 



+ V 2 n (k,K) (21) 



Note that for propagating waves, r\ n is real, i.e. if n > 0. If 
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uo n (tj = uo n (X) , (22) 

then obviously 

U n ± (kj = (kj ± r)n (X ^) J ( 2 3) 

which means that there is a gap, Au n , in the frequency spectrum 

Au n (k, k^j = u n+ (k, k u ^j - u n _ (k, k^j = 2r] n (k, k^ (24) 

A similar gap exists in the energy spectrum of electrons in a periodic potential 
[10] , with the following important difference: the width of the gap in the frequency 
spectrum of photons, Eq. (24), depends on the wave vector k. 

Equation (22) defines the surface that is the boundary of the Brillouin zone of 
the photon. For an isotropic medium, the frequency u n depends on the modulus 
of the wave vector, and the solution of Eq. (22) is k 2 = k 2 . Hence, the equation 
for the Brillouin zone boundary takes the form 

which is the equation of a plane in /c-space. Eq. (25) coincides with the well- 
known Bragg condition. 

Much more unusual is the photonic spectrum when two roots of Eq. (13) that 

belong to different branches (say, u n (jtj and u q (k^j) are close to each other, 

i.e. uj n {k^ ~ u q (k u ^j . Then, the solution of the dispersion equation Eq. (13) is 

given by Eq. (21) with u n (kuj replaced by u q (kjj and rj n (k, kjj substituted 
by 







0) 
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/i (k,u n f 2 (k u ,u n 




- U n 




^1 {k^ - LU q 



(26) 

To exhibit the boundary of the Brillouin zone, we consider small values of the 
wave vector, so that 



0Ji I k ) = 



Pi = 



d(k 2 ) 



i = n,q 



(27) 
(28) 
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The equality u n (jtj = 

— * 

surface in k— space: 



with 



uj q ( hi 



yields the following equation for the desired 



(k~k R y =r\ 

Z Pi 



(29) 



02 -Pi 



k v 



R 2 



P1P2 



W2-P1) 



kl + 



2'V 



^(O)-^i(O) 

P2-P1 



Equation (29) is the equations of a sphere with the center at k,R and having radius 
R. The spherical shape of the Brillouin zone originates from the interaction of 
different types of waves (n and q), and has no precedent in the theory of electrons 
in periodic lattices. In the absence of interactions, Eq. (29) loses its physical 
meaning, and R 2 becomes negative 

It is interesting to point out another distinctive feature of photonic periodic 
structures. In contrast to electrons, in electrodynamics the radiation (or scatter- 
ing) problem exists along with the eigenvalues problem considered above. In this 
case, the frequency, u, is an external parameter, and the dependence of the wave 
vector on the frequency, k (ui) , is to be found. This dependence can be obtained 
from the solution of the eigenvalues problem, u( k ). From Eqs. (16) and (23) it 
follows that 



When 



5k 



k 



k (u>) = k Q (cu) + 5 k (to). 



5u> — v g 5 k , 



(30) 
(31) 



where ~v g = ?m is the group velocity and 9 is the angle between v g and ko. 



dk 



3 Periodic array of resonant particles 

In this section, we consider a photonic crystal composed of periodically arranged 
identical dielectric particles having refractive index n, and study the spectrum 
of excitations (eigenwaves) associated with a resonance in a single particle. Such 
a resonance exists when the particle size, d, is of order of the wavelength inside 
the particle, A = — , (M i resonances, for example). This system is an analog of 
the tight-binding model in the solid state theory [12], and presents an opposite 
limiting case to that considered in the previous section. Indeed, while in a medium 
with weak periodic modulation, the band gap is small compared to the conduction 
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band, the gaps in the spectrum of the photonic crystal studied in this section are 
much wider than the passing zones. 

In what follows, we consider the frequency range in which the wavelength of 
the radiation in the medium (having refractive index no) is larger than the typical 
size of the particles 

d « — (32) 

UqUJ 

In this case, the problem can be handled by means of the local perturbation 
method (LPM) [13], [14], [15], [16]which is based on the assumption that the 
(unknown) field is independent of the coordinates inside each particle. In contrast 
to the Born approximation, which considers only weak scattering, LPM leaves 
room for arbitrarily large amplitudes of scattered fields, and does not rule out 
the existence of resonances that take place when d ~ — . 
Equation (1) now takes the form 



T)= — CXT) x 7 



(33) 



where the spatial-time dispersion of the homogeneous host medium is incorpo- 
rated in H. The operator U is non-zero only inside the particles and depends 
on their shape, position, and dielectric properties. For the sake of simplicity, we 
shall assume that the dielectric particles possess only time dispersion. 

Since U in Eq. (33) is a periodic function of r, the solution of Eq. (33) 
can be written in the form of Eq. (2). If the inequality (32) holds, the LPM 
approximation is applicable, implying that 

U(r- r p , to) E (r) e** &U(r- f p , to) E (r n ) e** ? * (34) 

where E (f) is a periodic function. If we assume that one of the particles is 
located at the origin, r — 0, then E (f p ) = E (0) . Substitution of Eq. (34) into 
Eq. (33) yields: 

o OO 

H( — ,cu)E(r) + E(0) ^ U (r - f p , u) = 0. (35) 

p=— oo 

Equation (35) can be rewritten in integral form, 

E(r) + E{0)Y^ I G{r-r)U(r-r p ,u)dr = 0, (36) 
v J 

where the Green function of the homogeneous (U = 0) host medium, G, is defined 
by the equation 

H(-?-,uj)G(r,r) = 5(r-r) (37) 
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Changing Eq. (37) to spatial Fourier transforms, we obtain 



1 r e ik (f-r) _^ 

G(f,r) = , / — =f d k (38) 

In an isotropic medium, H( k ,oj) — H(k 2 ,u), the integral in Eq. (38) can be 
calculated explicitly: 

G(f-r}= - - _ Y - 2 - , (39) 
where are roots of the equation 

H(-(l w)=0, (40) 

and ff= ^ 

By substituting r = in Eq. (36) we obtain the following dispersion equation 
for co (k) 

1 + / G(r)U(f,Lu)dr* + / U (r, uj)d?^G (r n ) e l ^ n = (41) 

It is apparent that the second term in Eq. (41) corresponds to the resonance 
frequencies of an isolated single particle, while the third term gives the corrections 
due to the multiple scattering (interactions). This term is small compared to the 
second one, therefore, we seek the solutions of Eq. (41) of the form 

u p = uj 0p + 5u p , (42) 

where u>o p are the solutions of the "zero-order" (without interactions) equation 

1 + J G(r,u)U(r,u)dr = 0, (43) 

and Su p « uop. 

Since the summands in Eq. (39) decrease exponentially with increasing \ f n \, 
we may keep only the term with the maximal d (denote it as ( m ). Then, for Su p , 
we find 

r f JU(r)df ^ exp(-C m \f n \) ^ 1 

P ~\~ £H\-C 2 J I G (r, u) U (r, u) dr ^ K\ * J ^ 

For the same reason, in the sum in Eq. (44), we take into account only the 
terms corresponding to the interaction with nearest neighbors. For a cubic crystal 
having side b, this yields 
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Su p = — t4(cos k x b + cos k y b + cos k z b), (45) 

A = JU(r) dfe~^m\b 

4nb£H'{-C 2 m )fG(?,u)U(?,u)dr 

The width of the conduction band is equal to 6 A « uj Po . To carry out the 
integration in Eq. (44), we consider the example of a homogeneous medium with 
time dispersion. Then the operator H in (35) is the Helmholtz operator 

H = A + ^e(u), (47) 
er 

where e(uS) is the dielectric permittivity of the host medium, in which the periodic 
grating of small particles (local perturbations) is located. For defmiteness, we 
assume that the dielectric constants of both host medium and particles are of the 
Drude form: 

£o M = 1 + ^-5; (48) 
£(") = !+ , 3 2 - (49) 

K) -u 2 

The dispersion equation Eq. (43) has two solutions: 



W l,2 - ^ ± 



n? + K) 2 + ui 



u?. K) 2 (50) 



where 

UJ 



2 47TC 2 



5 

is the frequency of the "geometrical" resonance related to the finite size, a, of 
the local perturbation (Debye resonance of the single particle); S = J f ^ dr ~a 2 . 
Note that the number of solutions depends on the explicit form of the functions 
£q(u) and s(oj), and can be larger than two. 

Multiple scattering between neighboring particles leads to the broadening of 
the discrete levels u 12 into passing zones of finite width. Frequencies ui 12 = 
^01,02 + ^1,2 correspond to a peculiar kind of polaritons that arise in photonic 
crystals with time dispersion due to the interaction of the electromagnetic field 
(iv), eigenoscillations of the dielectric medium (uq), and Debye resonance (uj s ). 

For the width of the interaction- induced passing zones, Awi ]2 , Eq. (45) yelds 

Awi,2 » 6 ^-^i, (51) 
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V= / f(r)df. 



(52) 



Terms independent on k (small shifts of iv l2 ) are omitted in Eq. (51). Recall 
that for the exponent in the Green function to be real (stability condition), it is 
nesessary that 

1 - — < 0. (53) 

^1,2 

One can see from (51) that the zone width is maximal at uii = Q and is equal 
Aojj = jUi. For —t — increasing, the width decreases exponentially, and at 

" u 1,2 

^ 2 >> 1, the zone structure disappears. When the spacing between levels, 

w 1,2 

\uji — uj 2 \ , is of order of u s , the width of the passing zone is - >> 1 times smaller 
than the size of the band gap. 

To conclude, time-spatial dispersion causes important differences in the spec- 
trum of photons in periodic dielectric systems from that of electrons in solid 
state crystals. It gives rise to different types of waves, changes the number and 
structure of passing zones and band gaps, and dramatically deforms the Brillouin 
zones. 
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